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rr 1 1 ABSTRACT 

o 

Context. In recent years, many low-mass (< 0.45M o ) white dwarf stars expected to harbor He cores have been detected in the field 
of the Milky Way and in several galactic globular and open clusters. Until recently, no objects of this kind showed pulsations. This 
situation has changed recently with the exciting discovery of SDSS J184037. 78+642312.3, the first pulsating low-mass white dwarf 
. star. 

Aims. Motivated by this extremely important finding, and in view of the very valuable asteroseismological potential of these objects, 
we present here a detailed pulsational study applied to low-mass He-core white dwarfs with masses ranging from 0.17 to 0.46M o , 
based on full evolutionary models representative of these objects. This study is aimed to provide a theoretical basis from which to 
interpret future observations of variable low-mass white dwarfs. 

Methods. The background stellar models on which our pulsational analysis was carried out were derived by taking into account the 
complete evolutionary history of the progenitor stars, with special emphasis on the diffusion processes acting during the white dwarf 
cooling phase. We computed nonradial g-modes to assess the dependence of the pulsational properties of these objects with stellar 
parameters such as the stellar mass and the effective temperature, and also with element diffusion processes. We also performed a 
g- and p-mode pulsational stability analysis on our models and found well-defined blue edges of the instability domain, where these 
stars should start to exhibit pulsations. 

Results. We found substantial differences in the seismic properties of white dwarfs with M, > 0.20M o and the extremely low-mass 
(ELM) white dwarfs (M, < 0.20M o ). Specifically, g-mode pulsation modes in ELM white dwarfs mainly probe the core regions and 
are not dramatically affected by mode-trapping effects by the He/H interface, whereas the opposite is true for more massive He-core 
white dwarfs. We found that element diffusion processes substantially affects the shape of the He/H chemical transition region, leading 
to non-negligible changes in the period spectrum of low-mass white dwarfs, in particular in the range of stellar masses characteristic 
of ELM objects. Finally, our stability analysis successfully predicts the pulsations of the only known variable low-mass white dwarf 
(SDSS J184037. 78+6423 12.3) at the right effective temperature, stellar mass and range of periods. 

Conclusions. Our computations predict both g- and p-mode pulsational instabilities in a significant number of known low-mass and 
ELM white dwarfs. It is worth observing these stars in order to discover if they pulsate. 
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O ! 1. Introduction Brown et al. 2010, 2012; Kilic et al. 2011, 2012); they are re- 

ferred to as extremely low-mass (ELM) white dwarfs. 

^ — 1 ' White dwarf stars are the most common evolutionary fate of 

^.1 low- and intermediate-mass (M, < 11M ; Siess 2007) stars. The low-mass white dwarf population is probably produced 

• rH ■ Indeed, the vast majority — more than 95% — of all stars will by strong mass-loss episodes at the red giant branch (RGB) 

1 die as white dwarfs. As such, these old and compact stellar rem- phase before the He-flash onset. As such, these white dwarfs are 

H . nants provide a wealth of information about the evolution of expected to harbor He cores, in contrast to average mass white 

. . ■ stars, star formation, and the age of a variety of stellar popula- dwarfs, which all likely contain C/O cores. For solar metallic - 

tions, such as our Galaxy and open and globular clusters (Winget ity progenitors (Z * 0.01 - 0.02), mass-loss episodes must oc- 

& Kepler 2008; Garcfa-Berro et al. 2010; Althaus et al. 2010). cur in binary systems through mass-transfer, since single star 

Among the different flavors of white dwarfs, the most common evolution is not able to predict the formation of these stars in 

is the spectral class of DA white dwarfs, characterized by H a Hubble time. This evolutionary scenario is confirmed by the 

rich envelopes, that comprises around 80% of all known white fact that most of low-mass white dwarfs are found in binary 

dwarfs. An important property of the white dwarf population systems (e.g., Marsh et al. 1995), and usually as companions 

is their mass distribution. For DA white dwarfs, the mass dis- to millisecond pulsars (van Kerkwijk et al. 2005). In particular, 

tribution peaks at w O.59M , and exhibits also high-mass and binary evolution is the most likely origin for ELM white dwarfs 

low-mass components (Kepler et al. 2007; Tremblay et al. 201 1; (Marsh et al. 1995). On the other hand, for high-metallicity pro- 

Kleinman et al. 2012). The population of low-mass white dwarfs genitors (Z « 0.03 - 0.05), the He-core flash can probably be 

has masses lower than 0.45 M and peaks at « O.39M . Recently, avoided by mass losses due to strong stellar winds on the RGB. 

a large number of low-mass white dwarfs with masses below This can be the origin of the isolated low-mass He-core white 

0.20 - O.25M has been discovered (Kawka & Vennes 2009; dwarfs (M, m O.4M ) found in the high-metallicity open cluster 
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NGC 6791 (Hansen 2005; Kalirai et al. 2007; Garcfa-Berro et al. 
2010). 

The evolution of low-mass white dwarfs is strongly depen- 
dent on their stellar mass and the occurrence of element diffu- 
sion processes. Althaus et al. (2001) and Panei et al. (2007) have 
demonstrated that element diffusion leads to a dichotomy re- 
garding the thickness of the H envelope, which translates into 
a dichotomy in the age of low-mass He-core white dwarfs. 
Specifically, for stars with M, > 0.18 - O.2OM , the white 
dwarf progenitor experiences multiple diffusion-induced ther- 
monuclear flashes that burn most of the H content of the en- 
velope, and as a result, the remnant enters the final cooling track 
with a very thin H envelope. In this way, the star is unable to 
sustain stable nuclear burning while cools and the evolution- 
ary timescale is rather short (« 10 7 yrfl On the other hand, if 
M„ < O.2OM , the white dwarf progenitor does not experience H 
flashes at all, and the remnant enters its terminal cooling branch 
with a thick H envelope. This is thick enough for residual H nu- 
clear burning to become the main energy source, that ultimately 
slows down the evolution, in which case the cooling timescale 
is of the order « 10 9 yrs. The age dichotomy has been also sug- 
gested by observations of those low-mass He-core white dwarfs 
that are companions to millisecond pulsars (Bassa et al. 2003; 
Bassa 2006). 

The internal structure of white dwarfs can be disentangled by 
means of asteroseismology (Winget & Kepler 2008; Fontaine & 
Brassard 2008; Althaus et al. 2010). White dwarf asteroseismol- 
ogy allows us to place constraints on the stellar mass, the thick- 
ness of the compositional layers, and the core chemical compo- 
sition, among other relevant properties. In connection with the 
core composition of white dwarfs, some theoretical work in the 
past has explored the differences that should be expected in the 
pulsation properties of low-mass white dwarfs harboring cores 
made either of C/O or He (Althaus et al. 2004), and high mass 
white dwarfs with cores made either of C/O or O/Ne (Corsico 
et al. 2004). More recently, Castanheira & Kepler (2008, 2009), 
and Romero et al. (2012) (see also Romero 2012) have carried 
out detailed asteroseismological studies on a large number of 
pulsating DA white dwarfs (DAV or ZZ Ceti variables), pro- 
viding valuable information about the range of the H envelope 
thicknesses expected for the class of DA white dwarfs. 

In order to apply the principles of asteroseismology on 
low-mass and ELM white dwarfs, it is necessary to find such 
stars (either in the Galactic field or in stellar clusters) under- 
going pulsations, something that until very recently had not 
been possible. In this regard, it is worth mentioning that ex- 
haustive searches reported by Steinfadt et al. (2012) have given 
null results, despite their theoretical predictions (Steinfadt et al. 
2010, hereinafter SEA 10) suggesting pulsational instabilities at 
least in stars with M„ < O.2OM (ELM white dwarfs). The 
situation has improved drastically with the exciting discovery 
of SDSS J184037.78+642312.3, the first pulsating ELM white 
dwarf (Hermes et al. 2012). This star (r eff = 9100 ± 170 K, 
logg = 6.22 + 0.06) exhibits multiperiodic photometric varia- 
tions with a dominant period at * 4698 s, much longer than the 
longest periodicities detected in any ZZ Ceti star (100 - 1200 
s). The discovery of this pulsating ELM white dwarf opens the 



1 Note that, if element diffusion processes are not included, theoret- 
ical computations also predict the occurrence of H-shell flashes before 
the terminal cooling branch is reached, but the H envelopes remain 
thick, with substantial H burning and large white dwarf cooling ages 
(Driebe et al. 1998; Sarna et al. 2000). 



possibility of sounding the interiors of low-mass white dwarfs 
employing the tools of asteroseismology. 



Needless to say, accurate and realistic stellar models rep- 
resentative of low-mass white dwarfs are needed to correctly 
interpret the present and future observations of pulsations in 
these objects. In particular, the complete evolutionary history of 
the progenitor stars must be taken fully into account in order 
to assess the correct thermo-mechanical structure of the white 
dwarf and the thickness of the H envelope remaining after multi- 
ple diffusion-induced CNO flashes occurring before the terminal 
cooling branch is reached, which strongly determines the cool- 
ing timescales. Also, element diffusion processes must be con- 
sidered in order to consistently account for the evolving shape 
of the internal chemical profiles (and, in particular, the chem- 
ical transition regions) during the white dwarf cooling phase. 
Finally, stable H burning, which is particularly relevant in the 
case of ELM white dwarfs (and can play an important role in the 
excitation of pulsations), must be taken into account as well. 



Motivated by the asteroseismic potential of pulsating low- 
mass white dwarfs, and stimulated by the discovery of the 
first variable ELM white dwarf, we have started in the La 
Plata Observatory a theoretical study of the pulsation proper- 
ties of low-mass, He-core white dwarfs with masses in the range 
0.17 - O.46M . Our interest in these stars is strengthened by the 
availability of detailed and fully evolutionary stellar models ex- 
tracted from the sequences of low-mass He-core white dwarfs 
presented by Althaus et al. (2009). These models were derived 
by considering the evolutionary history of progenitor stars with 
high metallicity, and taking into account a self-consistent, time- 
dependent treatment of the gravitational settling and chemical 
diffusion, as well as of residual nuclear burning. In this study, 
we explore the adiabatic pulsation properties of these models, 
that is, the expected range of periods and period spacings, the 
propagation characteristics of the pulsation modes, the regions 
of period formation, as well as the dependence on the effective 
temperature and stellar mass. For the first time, we assess the 
pulsation properties of He-core white dwarfs with masses in the 
range 0.20-0.45 M . In particular, we highlight the expected dif- 
ferences in the seismic properties of objects with M* > O.2OM 
and the ELM white dwarfs (M* < O.2OM ). Since the most im- 
portant work in the literature on pulsations of low-mass white 
dwarfs is the study of SEA 10, we shall invoke it repeatedly to 
compare their results with our own findings, in particular re- 
garding the seismic properties of ELM white dwarfs. We also 
explore the role of time-dependent element diffusion in ELM 
models on their pulsational properties, an aspect not discussed 
by SEA 10. In addition, we compute g- and p-mode blue edges 
of the instability domain of these stars through a nonadiabatic 
stability analysis of our sequences, and in particular, we try to 
determine whether our computations are able to predict the pul- 
sations exhibited by SDSS J184037.78+642312.3. This paper is 
organized as follows. In Sect. [2] we briefly describe our numeri- 
cal tools and the main ingredients of the evolutionary sequences 
we use to assess the pulsation properties of low-mass He-core 
white dwarfs. In Sect.[3]we present in detail our adiabatic pulsa- 
tion results, while in Sect.|4]we describe our nonadiabatic results 
and the case of SDSS J184037.78+642312.3. Finally, in Sect. [5] 
we summarize the main findings of the paper. 
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Fig. 1. The T e g - logg diagram showing our low-mass He-core white dwarf evolutionary tracks (black thin lines). Bold numbers 
correspond to the stellar mass of each sequence. The different symbols with error bars represent the location of the low-mass white 
dwarfs known to date. We include field objects as well as stars found in the open cluster NGC 6791 to have a general overview 
of the observational status of these stars. The location of the first known pulsating ELM white dwarf (T e ff = 9100 ± 170 K, 
log g = 6.22 + 0.06; Hermes et al. 2012) is displayed with a small circle (magenta). As reference, we also depict the C/O-core white 
dwarf evolutionary tracks (gray thick lines) for masses between 0.525 and O.935M (Romero et al. 2012) and the location of the 
known ZZ Ceti (DAV) stars. 



2. Computational tools and evolutionary sequences 

2.1. Evolutionary code and input physics 

The evolutionary models employed in our pulsational analysis 
were generated with the LPCODE evolutionary code, which pro- 
duces complete and detailed white dwarfs models incorporat- 
ing very updated physical ingredients. While detailed informa- 
tion about LPCODE can be found in Althaus et al. (2005, 2009) 
and references therein, we highlight here only those ingredients 
which are important for our analysis of low-mass, He-core white 
dwarf stars: 

(a) The standard mixing length theory (MLT) for convection 
(with the free parameter a = 1.6) has been adopted. Had we 
adopted a different convective efficiency, then the adiabatic 
pulsation results presented in this work (Sect. [3} should not 
be affected, although the predicted blue edge of the ELM and 
low-mass He-core white dwarfs instability domain (Sect. |4]i 
should appreciably change. 

(b) A supersolar metallicity for the progenitor stars has been 
considered: Z = 0.03. 

(c) Radiative opacities for arbitrary metallicity in the range from 
to 0.1 are from the OPAL project (Iglesias & Rogers 1996) 



supplemented at low temperatures with the molecular opaci- 
ties of Alexander & Ferguson (1994). 

(d) Neutrino emission rates for pair, photo, and bremsstrahlung 
processes were taken from Itoh et al. (1996), and for plasma 
processes we included the treatment of Haft et al. (1994). 

(e) Conductive opacities are those of Cassisi et al. (2007). 

(f) For the white dwarf regime we employed an updated version 
of the Magni & Mazzitelli (1979) equation of state. 

(g) The nuclear network takes into account 16 elements and 34 
nuclear reactions for pp chains, CNO bi-cycle, He burning, 
and C ignition. 

(h) Time-dependent diffusion due to gravitational settling and 
chemical and thermal diffusion of nuclear species was taken 
into account following the multicomponent gas treatment of 
Burgers (1969). 

(i) Abundance changes were computed according to element 
diffusion, nuclear reactions, and convective mixing. This de- 
tailed treatment of abundance changes by different processes 
during the white dwarf regime constitutes a key aspect in the 
evaluation of the importance of residual nuclear burning for 
the cooling of low-mass white dwarfs. 
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2.2. Pulsation codes 

In this work, we have carried out a detailed adiabatic g-mode 
pulsation study aimed at exploring the seismic properties of 
low-mass He-core white dwarf models. In addition, we have 
performed a g- and p-mode stability analysis of our model se- 
quences. The pulsation computations of our adiabatic survey 
were performed with the pulsation code described in detail in 
Corsico & Althaus (2006), which is coupled to the LPCODE 
evolutionary code. The pulsation code is based on a general 
Newton-Raphson technique that solves the full fourth-order set 
of equations and boundary conditions governing linear, adia- 
batic, nonradial stellar pulsations following the dimensionless 
formulation of Dziembowski (1971) (see Unno et al. 1989). The 
prescription we follow to assess the run of the Brunt-Vaisala fre- 
quency (AO for a degenerate environment typical of the deep in- 
terior of a white dwarf is the so-called "Ledoux Modified" treat- 
ment (Tassoul et al. 1990). 

As for the stability analysis, we have employed the finite- 
difference nonadiabatic pulsation code described in detail in 
Corsico et al. (2006). The code solves the full sixth-order com- 
plex system of linearized equations and boundary conditions as 
given by Unno et al. (1989). Typically, our adiabatic and non- 
adiabatic periods differ by less than « 2%. The caveat of our 
analysis is that the nonadiabatic computations rely on the frozen- 
convection approximation, in which the perturbation of the con- 
vective flux is neglected. While this approximation is known to 
give unrealistic locations of the g-mode red edge of instability, 
it leads to satisfactory predictions for the location of the blue 
edge of the ZZ Ceti (DAV) instability strip (see, e.g., Brassard 
& Fontaine 1999 and van Grootel et al. 2012) and also for the 
V777 Her (DBV) instability strip (see, for instance, Beauchamp 
et al. 1999 and Corsico et al. 2009). 



2.3. Evolutionary sequences 

To derive starting configurations for the He-core cooling se- 
quences consistent with the evolutionary history of the progeni- 
tor star, Althaus et al. (2009) simply removed mass from a 1M 
model at the appropriate stages of its evolution (Iben & Tutukov 
1986; Driebe et al. 1998). Other details about the procedure to 
obtain the initial models are provided in Althaus et al. (2009). 
The resulting final stellar masses (M*/M ) are listed in Table Q] 
with the total amount of H contained in the envelope (Mh/M*), 
and the coordinate - log(g) (being q — 1 - M r /M t ) of the loca- 
tion of the He/H chemical interface (where Xn e = Xn — 0.5), for 
models at r e ff » 10 000 K. The last column shows the time spent 
by the stars to cool from T eff ~ 14 000 K to * 8000 K. 

We mention again that a metallicity of Z = 0.03 for the 
progenitor stars has been adopted. This value is appropriate for 
overmetallic environments such as the open cluster NGC 6791 
(Althaus et al. 2009). Had we adopted a solar metallicity value 
(Z m 0.01 - 0.02) instead, then all our sequences would be 
characterized by H envelopes thicker than those considered here 
(second column of Table Q]), but it would not qualitatively affect 
the results presented in this paper. 

The sequence of M* = 0. 17M does not correspond to the set 
of sequences presented by Althaus et al. (2009); it was specif- 
ically computed for the present study. The models of this se- 
quence are representative of ELM white dwarfs. The evolution- 
ary history of this sequence is quite different as compared with 
the remaining sequences, because its progenitor star does not ex- 
perience any CNO flash. As a result, the remaining H envelope is 
markedly thicker than for the other sequences. The very distinct 



Table 1. Selected properties of our He-core white dwarf se- 
quences (metallicity of the progenitor star: Z = 0.03) at T e g ~ 
10 000 K: the stellar mass, the mass of H in the outer envelope, 
the location of the He/H chemical interface, and the time it takes 
the star to cool from T eS * 14 000 K to * 8000 K. 
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Fig. 2. The asymptotic period spacing in terms of the effective 
temperature for all of our low-mass He-core evolutionary se- 
quences. 



origin of the 0.1 7M Q sequence is also emphatically evidenced 
by the time the star spends to cool from w 14000 to * 8000 K, 
which is about 3.7 times longer than for the case of the 0. 198M 
sequence. 

The long time that ELM white dwarfs should spend at these 
effective temperatures, by virtue of vigorous stable H burning 
via the pp chain, motivated Steinfadt et al. (2012) to search for 
pulsating objects with masses lower than w O.2OM . We note 
from TableQ] however, that models with m 0.40-0.45 M should 
also have comparable evolutionary timescales. This implies that 
there should be a good chance of finding pulsating objects with 
these masses, if they exist, and current searches of variable low- 
mass white dwarfs should not be restricted only to ELM white 
dwarfs. 

In Fig. [1] we depict our evolutionary tracks on the plane 
r e (j - log g (black thin lines). In order to have an overview of the 
observational status, we have included the location of the low- 
mass white dwarfs known to date. They include field stars that 
are members of binary systems (Kilic et al. 2007, 201 1; Brown 
et al. 2012; Steinfadt et al. 2012, Silvotti et al. 2012) and sin- 
gle stars residing in the open cluster NGC 6791 (Kalirai et al. 
2007). As reference, we include also the location of the known 
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Fig. 3. The internal chemical profiles of He and H (upper panel) 
and the propagation diagram — the run of the logarithm of the 
squared critical frequencies — (lower panel) corresponding to an 
ELM white dwarf model of M. = O.17M and T eS * 9 000 K. 
Dots connected with thin dotted lines correspond to the spatial 
location of the nodes of the radial eigenfunction of low-order 
dipole {{ — 1) g-modes. 



ZZ Ceti stars and the evolutionary tracks corresponding to C/O- 
core white dwarfs (gray thick lines). 



3. Adiabatic pulsation properties 

3. 1 . Asymptotic period spacing 

For g-modes with high radial order k (long periods), the separa- 
tion of consecutive periods (\Ak\ = 1) becomes nearly constant at 
a value given by the asymptotic theory of nonradial stellar pul- 
sations. Specifically, the asymptotic period spacing (Tassoul et 
al. 1990) is given by: 



An? = n / 1), 
where 



r N - dr 



The squared Brunt-Vaisala frequency (AO is computed as: 

P X P 

where the compressibilities are defined as 
(d\nP\ (d\nP\ 



X P 



XT 



\dlnp) T{Xi) \dlnT) plXi} 
The Ledoux term B is computed as (Tassoul et al. 1990): 
d\nXi 



(1) 



(2) 



(3) 



(4) 



B 



XT V 



Xx, 



dlnP 



(5) 
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Fig. 4. Same as in Fig. |3] but for a M, = 0.45 M Q white dwarf 
model. 



where: 



Xx, 



d\nP 
dlnXi 



(6) 



P,T,(X i+i ) 



The expression in Eq. (Q} is rigorously valid for chemically ho- 
mogeneous stars. In this equation (see also Eq. [2]), the depen- 
dence on the Brunt-Vaisala frequency is such that the asymp- 
totic period spacing is larger when the mass and/or temperature 
of the model is lower. This trend is clearly evidenced by Fig. 
[2] for M, > O.25M , in which we depict the evolution of the 
asymptotic period spacing for all the sequences considered in 
this work. The higher values of All? for lower M» comes from 
the dependence N oc g, where g is the local gravity (g oc M*/Rl). 
On the other hand, the higher values of All? for lower T e ff result 
from the dependence oc yjxj, with xt — > for increasing 
degeneracy (T — > 0). 

Note that, however, the increase of All? with T e ff is not 
monotonous for the sequences with M* < O.25M . In fact, for 
this range of masses, the period spacing experiences a local min- 
imum before resuming the final growth trend with decreasing ef- 
fective temperature according to the predictions of Eq. @. That 
minimum in All? is found at lower T e ff as we go to lower masses. 
This somewhat weird behavior is induced by the evolving shape 
of the He/H chemical interface. We postpone an explanation for 
this feature of All? to Sect. [331 

The strong dependence of the period spacing on M, as evi- 
denced by Fig. |2]could be used, in principle, to infer the stellar 
mass of pulsating low-mass white dwarfs, provided that enough 
consecutive pulsation periods were available from observations. 
Such a prospect could be severely complicated by the fact that 
the period spacing of pulsating white dwarfs also depends on the 
thickness of the outer envelope (Tassoul et al. 1990), where All? 
is larger for thinner envelopes. This is particularly true in the 
context of low-mass He-core white dwarfs, in which the models 
of ELM objects harbor H envelopes that are several times thicker 
than more massive models. Then, regarding the value of An?, 
and for a fixed T e ff, a model with a low mass and a thick H enve- 
lope could readily mimic a more massive model with a thinner 
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Fig. 5. The run of the logarithm of the absolute value of the 
eigenfunctions y\ and y2, and the density of kinetic energy 
dEfkin/dr for t — 1 modes with k = 1 (upper panel), Ic — 5 
(middle panel), and k — 10 (lower panel), corresponding to the 
template model with M, = 0.1 7M . The vertical dashed line and 
the associated shaded strip mark the location of the He/H chem- 
ical transition region, whereas the vertical dotted line indicates 
the bottom of the outer convective zone. 

envelope. We envisage that, if a rich spectrum of observed peri- 
ods were available, this ambiguity could be broken by including 
additional information of the mode trapping properties, which 
yield clues about the H envelope thickness (see Sect. 13 .3b - 

3.2. Chemical profiles, critical frequencies and 
eigenfunctions 

Our low-mass white dwarf models are made of a He core sur- 
rounded by a H outer envelope. In between, there is a smooth 
transition region shaped by the action of microscopic diffusion, 
which is self consistently accounted for in LPCODE. In the up- 
per panels of Figs. [3] and |4] we display the internal chemical 
profiles for He and H corresponding to two template models of 
AT, = O.17M and M» = O.45M , respectively, at r eff * 9000 K. 
The less massive model is representative of ELM white dwarfs, 
and the more massive model is representative of massive He- 
core white dwarfs (0.20 < M„/M Q < 0.45). It is worth empha- 
sizing two important differences between these template models. 
First, the ELM model has a H envelope that is « 35 times thicker 
than the massive model. As mentioned in the Introduction, this 
is the result of the very different evolutionary history of the pro- 
genitor stars. Second, the other notorious difference between the 
chemical structures shown in Figs. [3] and [4] is the width of the 
He/H transition region. Specifically, this interface is markedly 
wider for the ELM model than for the massive one; see Section 
3.5 for an explanation. 

The features in the chemical profiles leave strong signatures 
in the run of the squared critical frequencies, in particular in A^ 2 , 
clearly shown in the lower panels of Figs. [3]and|4l which de- 
pict the propagation diagrams (Cox 1980; Unno et al. 1989) 
corresponding to these models. The squared Lamb frequency 
is defined as L 2 = ((( + l)c 2 /r 2 , where c s is the local adi- 
abatic sound speed, g-modes propagate in the regions where 
<x 2 < N 2 ,L 2 r where cr is the oscillation frequency. Note the 
very different shape of N 2 for both models. To begin with, the 



1 1 1 1 1 1 




| 1 | 1 | 1 | 1 | 1 | 1 


1 1 








k=1~ 


v He core 

1 i 1 i 




H envelope 

I i I i I i I i I i I i 


I | ' 


■ i 1 i 1 i 

"\ _/ 


.+/ 
II 


i i i i i i i i i i i i 

\-7<r>r ' ' /^^^^^^^^^^^-----^ 

f )' V l09ly1 ' f 

I | •— loglyj 


I 1 
k=5 - 


r v jj 
i i i 
' | 1 | 1 ' 




— log(dE„. /dr) 
i i i i km i i 


1 | ' 


1 l 1 l 1 l 




— — 1 — — 1 — — 1 — — 1 — — 1 — — 1 — 
I i I i I i I i I i I i 


1 1 










i 1! i ! 




1 1 ' m, = 0.45 m o 


k=10- 


■ii! 1 1 , 1 




1 , 1 , 1 , 1 , 1 , 1 , 


I i ■ 



I 1 I . I I I I I i I I I I I I I I I I I I I I I 

123456789 10 11 



-log(1-m/M.) 

Fig. 6. Same as Fig. |5] but for the model with M„ = 0.45 M Q . 

squared Brunt-Vaisala frequency for the ELM model is glob- 
ally lower than for the case of the O.45M model. This has to 
do with the lower gravity that characterizes the former model, 
which translates into smaller values of N 2 . Ultimately, a lower 
Brunt-Vaisala frequency profile produces a pulsation spectrum 
with longer pulsation periods. 

Another different feature between the propagation diagrams 
of both models is related to the bump at the He/H transition re- 
gion. This is notoriously more narrow and pronounced for the 
massive model than for the ELM one. As can be seen in Fig. [4] 
for the 0.45 M model this feature visibly affects the distribution 
of the nodes of the radial eigenfunctions, which cluster at the 
precise location of the peak in A^ 2 . In contrast, the distribution of 
nodes (and so, the shape of the eigenfunctions) is largely unaf- 
fected by the chemical interface in the case of the ELM model 
(Fig. 0. As we shall see in Sect. I3.3I this results in a weaker 
mode trapping as compared with the more massive models. 

Finally, another important characteristic of the ELM model 
is that its squared Brunt-Vaisala frequency exhibits a local max- 
imum at -log(g) * 2, which coincides with the loci of the 
bump due to the He/H transition region. But from there, A^ 2 in- 
creases very slightly outwards, without reaching much larger 
values in the surface layers. This is in contrast with the case 
of the O.45M model, in which the squared Brunt-Vaisala fre- 
quency exhibits larger values at the outer layers, resembling the 
situation in the C/O core DAV white dwarf models (see Fig. 3 
of Romero et al. 2012). This quite different shape of the run 
of Af 2 has strong consequences on the propagation properties 
of eigenmodes. Specifically, and due to the particular shape of 
Af 2 (that is, larger in the core than in the envelope; see also 
Fig. 2 of SEA 10), the resonant cavity of the eigenmodes for the 
ELM model is circumscribed to the core regions (- log q < 2), 
whereas for the O.45M , the propagation region extends along 
the whole model. 

This can be appreciated in more detail in Figs. [5] and [6] in 
which we plot the eigenfunctions y\ — £ r /r and j2 = (o~ 2 lg) £/, 
(where and £/, are the radial and the horizontal displace- 
ments, respectively) and the density of oscillation kinetic energy 
dE^in/dr (see Appendix A of Corsico & Althaus 2006 for its def- 
inition) corresponding to modes with radial order k — 1 , 5 and 
10 for the two template models (r eff a 9 000 K). From Fig. it 
is apparent that most of the spatial oscillations of the modes in 
the ELM white dwarf model are restricted to the regions below 
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tion kinetic energy (lower panel) in terms of the (t — 1) periods 
for the model with M. = O.17M and r eff = 9 000 K. The red 
horizontal line in the upper panel corresponds to the asymptotic 
period spacing computed with Eq. ([1]). 



the He/H interface, i.e., most of nodes are located in the region 
with - \og(q) < 2. This plot should be compared with Fig. 2 
of SEA 10. In summary, g-modes in ELM white dwarfs probe 
mainly the stellar corq3 and so, they have an enormous astero- 
seismic potential, as it was first recognized by SEA 10. 

In contrast, g-modes in more massive models (Hig. [6j oscil- 
late all along the stellar structure, and not just at the core, al- 
though they are more selectively affected by mode-trapping ef- 
fects produced by the compositional gradient at the He/H inter- 
face. In this sense, low-mass He-core white dwarfs with M» > 
O.2OM behave qualitatively similar to their massive cousins, the 
C/O-core DAV white dwarf stars. 

3.3. Mode trapping 

The period spectrum of chemically homogeneous stellar mod- 
els is characterized by a constant period separation, given very 
accurately by Eq. ([TJ- However, current evolutionary calcula- 
tions and different pieces of observational evidence indicate that 
white dwarf stars have composition gradients in their interiors 
(Althaus et al. 2010). The presence of one or more narrow re- 
gions in which the abundances of nuclear species (and so, the 
average molecular weight ft) are spatially varying modifies the 
character of the resonant cavity in which modes should propa- 
gate as standing waves. Specifically, chemical interfaces act like 
reflecting boundaries that partially trap certain modes, forcing 
them to oscillate with larger amplitudes in specific regions and 
with smaller amplitudes outside of those regions. The require- 
ment for a mode to be trapped is that the wavelength of its radial 
eigenfunction matches the spatial separation between two inter- 
faces or between one interface and the stellar center or surface. 
This mechanical resonance, known as mode trapping, has been 
the subject of intense study in the context of stratified DA and 
DB white dwarf pulsations (Brassard et al. 1992, Bradley et al. 
1993, Corsico et al. 2002a). In the field of PG 1 159 stars, mode 
trapping has been extensively explored by Kawaler & Bradley 
(1994) and Corsico & Althaus (2006). 



2 However, higher-order g-modes in the ELM model are less concen- 
trated in the core than low-order ones, and probe a significant portion 
of the star. 
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Fig. 8. Same as Fig.[7] but for the model with M„ = 0.45 M Q . 

There are two clear signatures of mode trapping. One of 
them affects the distribution of the oscillation kinetic energy 
(£kin) values in terms of the radial order (or the periods) of 
the modes. Since £kin is proportional to the integral of the 
squared amplitude of eigenfunctions weighted by the density 
(see Appendix A of Corsico & Althaus 2006), modes propagat- 
ing in the high-density environment typical of the deep interior 
of the star will exhibit higher values than modes that are oscil- 
lating in the low-density, external regions. Thus, a local maxi- 
mum in Efto is usually associated to a mode partially confined to 
the core regions, and a local minima corresponds generally to a 
mode partially trapped in the envelope. 

The second, and more important from an observational point 
of view, signature of mode trapping is that the forward period 
spacing An^ (= Tlk+i - HO, when plotted in terms of the pul- 
sation period n^, exhibits strong departures from uniformity. It 
is the period difference between an observed mode and adjacent 
modes (Ak = ±1) that matters as an observational diagnostic 
of mode trapping, at variance with E^, whose value is very 
difficult to estimate only from observations. For stellar models 
characterized by a single chemical interface, like the ones we 
are considering here, local minima in ATlk usually correspond to 
modes trapped in the outer layers, whereas local maxima in An^ 
are associated to modes trapped in the core region. 

In the upper panels of Figs. Q and [8] we show the forward 
period spacing for the template models with M, - O.17M and 
M» = 0.45 M , respectively. Lower panels depict the distribution 
of kinetic energy for the same models. The period spacings for 
the 0.1 7M model reach the asymptotic value of 103.72 s for pe- 
riods longer (» 5500 s) than the longest period shown in the plot. 
This is at variance with what happens in the case of the O.45M 
model, for which AYl^ approaches to the asymptotic value of 
64.21 s at periods m 1800 s. The signatures of mode trapping 
are much more notorious in the case of the O.45M model than 
for the O.17M one. Indeed, mode trapping in the ELM template 
model is appreciable only in the An* values for n* < 1500 s, 
and the E^n distribution is quite smooth for the complete range 
of periods shown. The weakness of mode-trapping in this model 
was anticipated in the previous section by virtue of the very wide 
and smooth He/H chemical transition region. In contrast, for the 
O.45M model the trapping signatures are stronger in both the 
ATlk and E-^m distributions, although the kinetic energy exhibits 
a quite smooth distribution for periods longer than * 2000 s. 

The fact that the asymptotic period spacing of the ELM 
model is not very rapidly reached (in terms of the pulsation peri- 
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Fig. 9. Pulsation periods of t — 1 g-modes in terms of the stellar 
mass for T e ff = 9500 K. Periods increase with decreasing M«. 



ods) is not in contradiction with the weakness of mode trapping 
effects characterizing this model. This is because, even for chem- 
ically homogeneous models — which completely lack of mode 
trapping effects due to the absence of chemical gradients — the 
asymptotic regime is reached for modes with radial orders higher 
than a given limit (minimum) value of k. This is clearly shown 
in Fig. 6 of Corsico & Benvenuto (2002), which displays the sit- 
uation for a chemically homogeneous white dwarf model made 
of pure He. 

In summary, our results indicate that mode trapping by the 
He/H interface is important in low-mass He-core white dwarfs 
with masses M, > O.2OM , but it is not an important issue in 
ELM white dwarfs. This feature severely limits the seismologi- 
cal potential of these objects to constrain the thickness of the H 
envelope. 

3.4. Dependence of the pulsations with Af. and r eff 



As discussed in Sect. 13.11 the period spacing and the periods 
vary as the inverse of the Brunt-Vaisala frequency. This critical 
frequency, in turn, increases with larger stellar mass and with 
higher effective temperature. As a result, the pulsation periods 
are longer for smaller mass (lower gravity) and lower effective 
temperature (increasing degeneracy). The influence of the stel- 
lar mass on the pulsation periods is shown in Fig. [9] where we 
plot Hi for I = 1 in terms of M* for a fixed value of r e ff . Note 
that the increase of with decreasing M* is less pronounced 
for low-order modes than for higher-order ones. This is in agree- 
ment with the behavior typical of the ZZ Ceti star models (see 
Bradley 1996). In Fig.[10]we show the evolution of the pulsation 
periods with T e s for the case of models with M„ = 0.30M o . The 
lengthening of the periods as the effective temperature drops is 
evident, although the effect is less pronounced than the decrease 
with the stellar mass (as compared with Fig. [9]). 



The effect of varying M, on the period-spacing and kinetic- 
energy distributions for a fixed T e ff is displayed in Fig. [TT] in 
which we present the results for three white dwarf models with 
stellar masses of M, = 0.17,0.30 and O.45M and the same ef- 
fective temperature (r e jf w 9000 K). As expected, AYl^ strongly 
changes with M», as An a does. In fact, the average of the An^ 
values varies from w 60 s to as 100 s (« 60%) when the stel- 
lar mass changes from 0.45 M Q to O.17M . Also, substantial 
changes in the kinetic energy distribution are expected when we 
vary the stellar mass, as evidenced by the lower panel of the Fig. 

Finally, we examine the effect of changing T e s on An* and 
£kin with fixed M* in Fig. [12] where we show these quantities in 
terms of the pulsation periods for a model with M* = O.3OM 
and three different values of the effective temperature. The sen- 
sitivity of the period spacing and the kinetic energy on r e ff are 
negligible for short periods (II* < 500 s) but it is quite appre- 
ciable for longer periods. Specifically, the average of the An* 
values varies from » 70 s to « 85 s (« 20%) when the effective 
temperature changes from 12 000 K to 8000 K. Thus, the depen- 
dence of AH* with r e ff is substantially weaker than with M,. The 
kinetic energy values also are quite sensitive to changes in r e tf, 
at least for periods longer than m 500 s. For the hottest model de- 
picted in Fig.[l2](7 , efF = 12 000 K), the more energetic modes are 
those with k = 1, ■ • • , 5, because they oscillate deep in the star 
where the density is high. However, for the cooler models (9500 
and 8000 K), the kinetic energy of the high-order modes, which 
are characterized by eigenfunctions mainly concentrated in the 
outer layers, strongly increase. This is because the decrease in 
T e ff produces a deepening of the outer convection zone, and so, 
the high-order modes feel gradually the presence of convection 
as the white dwarf cools down. Since g-modes in a convection 
zone become evanescent, such modes are forced to gain larger 
amplitudes at regions somewhat below the base of the outer con- 
vection zone, where the density is larger. Since E^ n is propor- 
tional to the integral of the squared eigenfunctions, weighted 
by p, these modes oscillate with larger energies. However, low- 
order modes remain rather insensitive to the thickening of the 
outer convection zone. 



3.5. Effects of element diffusion 

Finally, we discuss the effects that evolving chemical profiles 
have on the pulsation properties of low-mass, He-core white 
dwarfs. The effect of element diffusion on DAV stars has been 
explored by Corsico et al. (2002b). Time-dependent element dif- 
fusion, which is fully taken into account in our computations, 
strongly modifies the shape of the He and H chemical profiles 
as the white dwarf cools, causing H to float to the surface and 
He to sink down. In particular, diffusion not only modifies the 
chemical composition of the outer layers, but also the shape of 
the He/H chemical transition region itself. This is clearly borne 
out from Fig. [l3]for the case of the O.17M sequence in the T e g 
interval (1 1 000 - 8000 K). For the model at r eff = 11 000 K, 
the H profile is characterized by a diffusion-modeled double- 
layered chemical structure, which consists in a pure H enve- 
lope atop an intermediate remnant shell rich in H and He (up- 
per panel). This structure still remains, although to a much less 
extent, in the model at r e (j = 9500 K (middle panel). Finally, 
at r e ff = 8000 K, the H profile has a single-layered chemical 
structure (lower panel). Element diffusion processes affect all 
the sequences considered in this paper, although the transition 
from a double-layered structure to a single-layered one occurs 
at different effective temperatures. The markedly lower surface 
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Fig. 10. Pulsation periods of £ — 1 g-modes in terms of the ef- 
fective temperature for M* = O.3OM . Periods increase with de- 
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Fig. 11. The dipole forward period spacing (upper panel) and 
the oscillation kinetic energy (lower panel) in terms of the pe- 
riods for models with M t = 0.17 M Q , M, = O.3OM , and 
M t = O.45M , all of them characterized by a 7/ eff = 9 000 K. 
The horizontal lines in the upper panel correspond to the asymp- 
totic period spacings computed with Eq. (Q]). 



gravity that characterizes the less massive model, results in a less 
impact of gravitational settling, and eventually in a wider chem- 
ical transition in the model with M„ = O.17M . Because of this 
fact, the sequences with larger masses reach the single-layered 
structure at higher effective temperatures. 

The changes in the shape of the He/H interface are translated 
into non-negligible changes in the profile of the Brunt-Vaisala' 
frequency, as can be appreciated in the plot. In fact, at high ef- 
fective temperatures, N is characterized by two bumps, which 
merge into a single one when the chemical profile at the He/H in- 
terface adopts a single-layered structure. The fact that the chem- 
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Fig. 12. Same as Fig. Q~T] but for the case of a model with 
M, = O.3OM and three different effective temperatures of 
r eff = 12 000, 9500 and 8000 K. 



ical profiles evolve with time in our models is at variance with 
the claim by SEA 10, who argue that in their ELM models, dif- 
fusive equilibrium is valid in the He/H transition region. In the 
light of our results, we conclude that this assumption is wrong. 

In order to quantify the impact of the diffusively evolving 
chemical profiles on the period spectrum of low-mass He-core 
white dwarfs, we focus on the same ELM sequence analyzed 
in Fig. Q2] Specifically, we have considered an additional se- 
quence of models for which we switch off chemical diffusion 
in the LPCODE evolutionary code from T s $ m 1 1 000 K down- 
wards. This allow us to isolate the effect that diffusion processes 
induce on the pulsation periods from those changes induced by 
cooling. We found substantial changes in the value of the peri- 
ods when diffusion is neglected, depending on the specific mode 
considered and the value of the effective temperature. For in- 
stance, for the k = 3 mode (n 3 * 460 s) at 7/ eff * 9000 K, a 
variation of 5% in the value of the period is expected. 

We conclude that time-dependent element diffusion does af- 
fect the pulsation spectrum of low-mass and ELM white dwarfs, 
and that it must be taken into account in future pulsational anal- 
ysis of ELM white dwarfs. 

We close this section with an explanation to the behavior of 
An^ in terms of the effective temperature for the sequences with 
M, < O.22M , shown in Fig. [2] (Sect. 13. U . As we have seen, 
element diffusion processes lead to the formation of a double- 
layered configuration for the H profile in our models. This con- 
figuration consists in two regions where the He and H abun- 
dances exhibit a steeped spatial variation, which translates into 
two local features in the Brunt-Vaisala frequency. As evolution 
proceeds, the outer interface shifts downwards, gradually getting 
closer to the internal interface. The merging of the two bumps in 
N tends to increase the integrand of Eq. (0, thus producing a de- 
crease in the asymptotic period spacing. Thus, for a given T e ff, 
reaches a local minimum (see Fig. |2}. Later, cooling again 
dominates and then An^ increases, as expected. We found that 
the effective temperature at which AFL^ reaches the minimum is 
higher for larger stellar masses. Thus, this minimum is not seen 
in Fig.[2]for the sequences with M„ > O.25M because in these 
cases that minimum lies out of the range of T e ff considered. 
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arithm of the squared Brunt-Vaisala frequency, for a model with 
M, = 0.1 7M at different effective temperatures, which are in- 
dicated in each panel. 




Fig. 14. A zoom of the r e ff-log g diagram shown in Fig.Q]show- 
ing the region where is located the only known (so far) pulsating 
ELM white dwarf, SDSS J184037.78+642312.3 (Hermes et al. 
2012) displayed with a red circle. We include the approximate 
I — 1 g-mode blue edge of the theoretical ELM instability do- 
main derived by SEA 10 (thick black solid line), and the t = 1 
g-mode blue edge computed in this work (blue solid line) for 
all the sequences considered, included that of M* = O.17M , 
and extended towards lower gravities (blue dashed line). We also 
show our I = 1 p-mode blue edge of instability (thin dotted line). 



4. Nonadiabatic computations and the case of 
SDSS J1 84037.78+64231 2.3 

Having thoroughly described in the previous sections the g- 
mode adiabatic seismic properties of our low-mass He-core 
white dwarf models, in what follows we will focus our atten- 
tion on the only (up to now) ELM white dwarf discovered to 
be apulsating star, SDSS J184037.78+642312.3.This DA white 
dwarf, which is the coolest (T e s = 9100 + 170) and the lowest- 
mass (M, » 0.1 7M ) known pulsating (1200 < n < 4500 s) 
object of this kind, was discovered by Hermes et al. (2012) in 
the context of a systematic search for variable He-core white 
dwarfs among the large number of ELM white dwarfs discov- 
ered through the ELM Survey (Brown et al. 2010, 2012; Kilic et 
al. 2011, 2012). The discovery of SDSS J184037.78+642312.3 
opens the possibility to use the tools of white dwarf asteroseis- 
mology to obtain valuable information about the internal struc- 
ture of low-mass white dwarf stars. 

Here, we do not attempt an asteroseismological fitting to 
this star, because the period data are still far from be definitive. 
Rather, we concentrate mainly on the question of whether our 
low-mass He-core white dwarf models are able to predict the 
pulsations exhibited by SDSS J184037.78+642312.3 at the right 
effective temperature, stellar mass and range of observed peri- 
ods. 

In order to investigate the plausibility of excitation of pul- 
sations in our models, we performed a linear stability analysis 
on our complete set of evolutionary sequences. We found that a 
dense spectrum of g-modes are excited by the k - y-mechanism 
acting in the H partial ionization zone for all the masses con- 
sidered, and that there exists a well-defined blue (hot) edge of 



instability of He-core white dwarfs, which is the low-mass ana- 
log to the blue edge of the ZZ Ceti instability strip. We also 
found numerous unstable p-modes; the corresponding blue edge 
is hotter than the g-mode blue edge and has a lower slope. In 
Fig. [14] we show a zoom of the Fig.Q]at the region where SDSS 
J184037.78+642312.3is located. We have included the location 
of the theoretical dipole (I — 1) blue edges of the instability 
domain of low-mass white dwarfs. The g-mode blue edge com- 
puted in this work is marked with a blue solid line, and it is 
extended towards lower gravities with a dashed blue line. The p- 
mode blue edge is displayed with a thin dotted line. Interestingly, 
SDSS J184037.78+642312.3 is slightly cooler than the derived 
g-mode blue edge, and so, our computations does predict long- 
period pulsations in this star. It should be kept in mind, however, 
that this result could change if a MLT prescription with less con- 
vective efficiency were assumed. Our stability analysis relies on 
the frozen-convection approximation that is known to give sat- 
isfactory predictions for the location of the g-mode blue edge of 
other pulsating white dwarfs (DAVs, see Brassard & Fontaine 
1999; DBVs, see Corsico et al. 2009). Recently, a detailed study 
by van Grootel et al. (2012) showed that the blue edge of DAV 
models with O.6M computed with the frozen convection ap- 
proximation does not dramatically differ from that obtained with 
a time-dependent convection treatment in the nonadiabatic com- 
putations. We note that also p-modes are predicted by our anal- 
ysis to be destabilized in SDSS J184037.78+642312.3, although 
no short-period luminosity variations characteristic of acoustic 
modes have been detected in this star so far. 

We emphasize that, besides SDSS J184037.78+642312.3, 
our nonadiabatic analysis predicts g- and /?-mode pulsational in- 
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Table 2. g-mode pulsation quantities of the ELM white 
dwarf model representative of the pulsating star SDSS 
J184037.78+642312.3, characterized by M, = O.17M , r eff = 
9099 K, logg = 6.28, and log(L/L ) = 0.0156. 
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Fig. 15. Upper panel: the instability domain (blue dots) on the 
T e s - II plane for ( — 1 g-modes corresponding to our set of 
models with M, - 0.17M o . Also shown are the periodicities 
measured in SDSS J184037.78+642312.3 (red points with error 
bars). Lower panel: same as upper panel, but for the p-mode 
instability domain (black dots). 



stabilities in the other stars at the right of the blue edges. We 
believe that it would be worthwhile to examine photometrically 
these stars in order to see whether they are pulsating or not. 

The thick black solid line in Fig. [T4l corresponds to the g- 
mode blue edge of the theoretical ELM white dwarf instability 
region derived by SEA 10 by using the instability criterion II < 
Stt T th (Brickhill 1991; Wu & Goldreich 1999) satisfied by i = 1 
and k — 1 modes, where t± is the thermal timescale at the base 
of the convection zone. This criterion has been shown to be accu- 
rate in the case of ZZ Ceti pulsators, which are characterized by 
significantly higher gravities, but SEA10 were forced to extrap- 
olate its predictions to the low-gravity regime of the ELM white 
dwarfs. The blue edge of SEA10 is about 350 K cooler than ours 
at logg » 6.2, and in addition, its slope is somewhat larger, most 
likely the result of their using a different parametrization of con- 
vective efficiency (a = 1.0, Steinfadt 2011). At variance with 
our results, the blue edge derived by SEA 10 fails to predict pul- 
sational instabilities in SDSS J184037.78+642312.3. We note 
that, in order to estimate the blue edge, SEA 10 implicitly assume 
that the k = 1, { = 1 mode must be unstable, something that is 
not guaranteed at the outset without performing detailed nonadi- 
abatic pulsation computations (see below). If the blue edge were 
derived by considering higher order modes (k > 1), then the dis- 
crepancy would be even larger because their blue edge should 
shift to lower effective temperatures. 

We explored the domain of unstable dipole modes in terms 
of the effective temperature for the sequence of M„ = O.17M . 
For this sequence, g-modes become unstable at T e s as 9200 K. 
Fig. Q3] show the instability domains of both g- and p-modes 
on the T e g - II diagram. In the case of g-modes, we found that 
unstable modes have radial orders k > 9 (II > 1100 s), where 
the most unstable modes for each value of r e ff have the longest 
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period. Modes with k - 6,1 and 8 are only marginally unsta- 
ble (their stability coefficients are extremely small) and have 
not been included in the plot. Finally, modes with k - 1 , ■ ■ ■ , 5 
are pulsationally stable. In the case of p-modes, we found that 
these modes destabilize at somewhat higher effective tempera- 
tures (r e ff « 9330 K). At variance with the case of g-modes, 
low-order p-modes with k = 1,2,3, ■ • 29 (109 > II > 7.5 s) 
simultaneously become unstable at the blue edge, being the most 
unstable ones (largest stability coefficients) those with the short - 
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est periods. For the unstable models, the asymptotic frequency 
spacing of p-modes, defined as (Unno et al. 1989): 

ranges from Av a = 4.67 mHz (at T eff « 9330 K) to Av a = 6.20 
mHz (at T eff * 8000 K) for our M» = O.17M sequence. 

We have completed a reanalysis of the periods detected in 
SDSS J184037.78+642312.3 using nearly 33.0 hr of photomet- 
ric observation from October 2011 to July 2012. Our observ- 
ing and reduction techniques are identical to those described in 
Hermes et al. (2012), but include more than 27.0 hr of additional 
observations. Aliasing of the data has made determining the pe- 
riods especially difficult, so we have made every attempt to be 
as conservative as possible in our error estimates. We include in 
Fig.[l5]the periods detected in SDSS J184037.78+642312.3 as 
red dots with their corresponding uncertainties. 

The highest-amplitude mode excited in the star occurs at 
4697.8 ±4.3 s (51 +5 mma); this periodicity appears in the 5.5 hr 
of October 201 1 discovery data, but was the second-highest alias 
of the 4445 s mode quoted by Hermes et al. (2012). Additionally, 
we see evidence of variability at the following periods: 4310 ± 
200 s (28 ±3 mma), 2309+60 s (25 + 3 mma), 2400± 120 s (21 ±3 
mma), 4890 + 270 s (2 1 + 3 mma), and 2094 ± 50 s ( 1 7 + 3 mma). 
These periodicities may be independently excited modes, but it 
is possible that the shorter-period modes are nonlinear combi- 
nation frequencies of the longer modes. Additionally, this ELM 
white dwarf has an unseen companion; spectroscopic observa- 
tions by Brown et al. (2012) find it to be in a 4.6 hr binary 
with a minimum O.64M companion. If the ELM white dwarf 
is synchronized (or nearly synchronized) with the orbital period, 
it is possible that the modes just short and long of the 4697.8 s 
mode may be rotationally split components rather than indepen- 
dent modes. 

Our g-mode stability computations are in good 
agreement with the range of periods observed in SDSS 
J184037.78+642312.3, as shown in Fig.[l5] The small discrep- 
ancy for the longest periods can easily be accommodated taking 
into account the uncertainties in T e ff, which are known to be 
severe because the uncertainties characterizing the atmosphere 
models of very low-mass white dwarfs (Steinfadt et al. 2012). 

Finally, we consider an ELM white dwarf model represen- 
tative of SDSS J184037.78+642312.3. From our suite of mod- 
els, we found a 0.1 7M white dwarf model with T e ff = 9099 
K, logg = 6.28, and log(L/L Q ) = 0.0156. The asymptotic pe- 
riod spacing of g-modes is An a = 103.47 s, and the asymptotic 
frequency spacing for p-modes is Av a = 4.87 mHz. In Table 
|2] we show the g-mode periods, period spacings, rates of pe- 
riod change and oscillation kinetic energies corresponding to this 
model. Note that the longest observed period (« 4700 s) should 
correspond in our model to a I — 1 g-mode with a very high 
radial order, k =s 44 - 45. By examining the radial eigenfunction 
and density of kinetic energy of these modes, we found that most 
(=s 64%) of their nodes lay within the He core, although a sig- 
nificant number of nodes (16) are located at regions beyond the 
He/H interface. In other words, the observed mode with period 
=s 4700 s would sample a significant portion of the interior of the 
star, besides the stellar core. The I — 1 asymptotic period spac- 
ing, of a; 103 s is more than two times longer than for a normal 
O.6M ZZ Ceti star at r eff « 12 000 K. So, in addition to the low 
gravities and effective temperatures, the pulsational signatures of 
ELM white dwarfs should make them easy to distinguish them 
from the well-known, C/O-core ZZ Ceti stars. 



5. Summary and conclusions 

We have presented in this paper the first comprehensive the- 
oretical study of the seismic properties of low-mass, He-core 
white dwarfs with masses in the range 0.17-0.46MeB We have 
employed fully evolutionary stellar structures representative of 
these stars, extracted from the sequences of low-mass He-core 
white dwarfs published by Althaus et al. (2009). These mod- 
els were derived by considering the evolutionary history of pro- 
genitor stars with supersolar metallicities, and accounting for a 
time-dependent treatment of the gravitational settling and chem- 
ical diffusion, as well as of residual nuclear burning. We have 
explored the adiabatic pulsation properties of these models, in- 
cluding the expected range of periods and period spacings, the 
propagation properties and mode trapping of pulsations, the re- 
gions of period formation, as well as the dependence on the ef- 
fective temperature and stellar mass. For the first time, we as- 
sessed the pulsation properties of He-core white dwarfs with 
masses in the range 0.20 - O.45M . In particular, we empha- 
size the expected differences in the seismic properties of objects 
with AT. > 0.20M o and the ELM white dwarfs (M. < O.2OM ). 
The pulsation properties of ELM white dwarfs have been ex- 
plored recently by SEA 10. Our work can be considered as com- 
plementary to that study. We have also explored the role of time- 
dependent element diffusion in ELM white dwarf models on 
their pulsational properties. Finally, we have computed g- and p- 
mode blue edges of the instability domain for these stars through 
a nonadiabatic stability analysis. In particular, we attempted to 
determine if our evolutionary/pulsation models are able to pre- 
dict the pulsations exhibited by SDSS J184037.78+642312.3 
both in the right effective temperature and mass, and also in the 
correct range of periodicities. 

We summarize our findings below: 

- Extremely low mass (ELM) white dwarfs (M, < O.2OM ) 
have a H envelope that is much thicker than for massive He- 
core white dwarfs (« 0.20 - O.45M ), due to the very differ- 
ent evolutionary history of the progenitor stars. In particu- 
lar, ELM white dwarfs did not experience diffusion-induced 
CNO flashes, thus they harbor thick H envelopes. 

- By virtue of the thicker H envelope characterizing ELM 
white dwarfs, they experience H burning via the pp chain, 
and consequently their evolution is extremely slow. This fea- 
ture makes these stars excellent candidates to become pulsat- 
ing objects (SEA10). However, He-core white dwarfs with 
M* « 0.40 - 0.45 M Q should have evolutionary timescales of 
the same order, also making them attractive targets for cur- 
rent searches of variable low-mass white dwarfs. 

- The thickness of the He/H transition region is markedly 
wider for the ELM white dwarfs than for massive He-core 
white dwarfs. This is due to the markedly lower surface grav- 
ity characterizing ELM white dwarfs, that results in a less 
impact of gravitational settling, and eventually in a wider 
chemical transition. 

- The Brunt-Vaisala' frequency for ELM white dwarfs is glob- 
ally lower than for the case of massive objects, due to the 
lower gravity characterizing ELM white dwarfs. A lower 
Brunt-Vaisala frequency profile leads to longer pulsation pe- 
riods. 



3 Detailed tabulations of the the stellar models, chemical profiles, and 
pulsation periods of g- and p-modes for different stellar masses and 
effective temperatures are available at our web site http://www.fcaglp 
unlp.edu. ar/evolgroup. 
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The bump of N 2 at the He/H transition region is notoriously 
more narrow and pronounced for the massive white dwarfs 
than for the ELM white dwarfs. This feature results in a 
weaker mode trapping in ELM white dwarfs, something that 
severely limits their seismological potential to constrain the 
thickness of the H envelope. 

As already noted by SEA 10, the Brunt- Vaisala frequency 
of the ELM white dwarfs is larger in the core than in the 
envelope. This is in contrast with the case of models with 
M > O.2OM , in which the Brunt-Vaisala frequency exhibits 
larger values at the outer layers, thus resembling the situa- 
tion encountered in ZZ Ceti stars. So, g-modes in ELM white 
dwarfs probe mainly the stellar core and have an enormous 
asteroseismic potential, as it was first recognized by SEA 10. 
Similarly to ZZ Ceti stars, the g-mode asymptotic period 
spacing in low-mass He-core white dwarfs is sensitive pri- 
marily to the stellar mass, and to a somewhat less extent, 
to the effective temperature. Specifically, An 3 , is longer for 
lower M, and T e ff. Also, there is a non-negligible depen- 
dence with the thickness of the H envelope, where An 3 , is 
longer for thinner H envelopes (small Mh). Typically, the 
asymptotic period spacing range from m 55 s for M t = 
0.45 M and T eS = 11 500 K, up to « 1 10 s for M« = O.17M 
and T e ff = 8000 K (see Fig. [2j. The evolution of the asymp- 
totic period spacing in terms of the effective temperature is 
markedly influenced by the evolving chemical profiles due 
to element diffusion. 

Similar to what happens with the asymptotic period spacing, 
the spectrum of g-mode pulsation periods themselves is very 
sensitive to the stellar mass, the effective temperature and the 
thickness of the H envelope. Again, the largest dependence 
is with the stellar mass. 

Time-dependent element diffusion does appreciably affect 
the g-mode pulsation spectrum of low-mass white dwarfs, 
in particular ELM white dwarfs. In this regard, it is expected 
that diffusion processes substantially modify the He/H chem- 
ical interface, and so the resulting period spectrum, for mod- 
els with stellar masses M» < O.22M . We emphasize that 
time-dependent element diffusion should be taken into ac- 
count in future pulsational analysis of these stars. 
We obtained for the first time the g-mode blue (hot) edge 
of the instability domain for He-core low-mass white dwarfs 
with 0.20 < M»/M < 0.45. On the other hand, our g-mode 
blue edge for ELM white dwarfs (M t < O.2OM ) is hotter 
than that found by SEA 10. 

We also found a p-mode blue edge of instability of low-mass 
He-core white dwarfs. So, according to our analysis, several 
known low-mass white dwarfs with masses below ss O.3OM 
should be pulsating stars showing short and long periods as- 
sociated to p- and g-modes, respectively (see Fig. IT4l>. 
Our stability analysis successfully predicts the pulsations 
oberved in the only known variable low-mass white dwarf, 
SDSS J184037.78+642312.3 (Hermes et al. 2012), at the 
right effective temperature, stellar mass and range of pul- 
sation periods. Our computations also predict the presence 
of short-period pulsations in this star that, however, have not 
been detected by now. 

We found a representative model of SDSS 
J184037.78+642312.3, with parameters AT. = 0.1 7M Q , 
7/ eff = 9099 K, logg = 6.28, log(L/L ) = 0.0156, 
An a = 103.47 s and Av a = 4.87 mHz. According to this 
model, the longest pulsation period exhibited by this star 
(=s 4700 s) would be associated to a t = 1 g-mode with a 



very high radial order, k — 44 — 45, which probably samples 
a significant portion of the star, apart from the core. 

The recent discovery of SDSS J184037.78+642312.3, the 
first pulsating low-mass white dwarf star, has opened a new op- 
portunity to sound the interiors of these kind of stars through 
asteroseismology. In the next years a significant number of these 
pulsating objects will likely be uncovered through systematic 
photometric searches like the ones performed by Steinfadt et 
al. (2012) and Hermes et al. (2012). Needless to say, in order 
to accurately decode the information embedded in the pulsation 
spectrum of these stars, it will be necessary to have at hand a 
large suite of detailed evolutionary/pulsational models of low- 
mass white dwarfs. This paper is aimed at fulfilling this require- 
ment. 

Acknowledgements. Part of this work was supported by AGENCIA through the 
Programa de Modernization Tecnologica BID 1728/OC-AR, and by the PIP 1 12- 
200801-00940 gTant from CONICET. This research has made use of NASAs 
Astrophysics Data System. 



References 

Alexander, D. R., & Ferguson, J. W. 1994, ApJ, 437, 879 

Althaus, L. G., Corsico, A. H., Isern, J., & Garcfa-Berro, E. 2010, A&AR, 18, 
471 

Althaus, L. G., Panei, J. A., Romero, A. D., et al. 2009, A&A, 502, 207 
Althaus, L. G., Serenelli, A. M., Panei. J. A., et al. 2005, A&A, 435, 631 
Althaus, L. G., Corsico, A. H., Gautschy, A., et al. 2004, MNRAS, 347, 125 
Althaus, L. G., Serenelli, A. M., & Benvenuto O. G., 2001, MNRAS, 323, 471 
Bassa, C. G., 2006, PhD Thesis, Universiteit Utrecht 

Bassa, C. G., van Kerkwijk, M. H., & Kulkarni, S. R., 2003, A&A, 403, 1067 
Beauchamp, A., Wesemael, F., Bergeron, P., et al. 1999, ApJ, 516, 887 
Bradley, P. A. 1996, ApJ, 468, 350 

Bradley, P. A., Winget, D. E., & Wood, M. A. 1993, ApJ, 406, 661 

Brassard, P., & Fontaine, G. 1999, Stellar Structure: Theory and Test of 

Connective Energy Transport, 173, 329 
Brassard, P., Fontaine, G., Wesemael, F., & Hansen, C. J. 1992, ApJS, 80, 369 
Brickhill, A. J. 1991, MNRAS, 251, 673 

Brown, W. R., Kilic, M., Allende Prieto, C, & Kenyon, S. J. 2012, ApJ, 744, 
142 

Brown, W. R., Kilic, M., Allende Prieto, C, & Kenyon, S. J. 2010, ApJ, 723, 
1072 

Burgers J. M., 1969, "Flow Equations for Composite Gases" (New York: 
Academic Press) 

Cassisi, S. Pothekin, A. Y., Pietrinferni, A., Catelan, M. & Salaris, M., 2007, 
ApJ. 661, 1094 

Castanheira, B. G., & Kepler, S. O. 2009, MNRAS, 396, 1709 

Castanheira, B. G., & Kepler, S. O. 2008, MNRAS, 385, 430 

Corsico, A. Ft., Althaus, L. G., Miller Bertolami, M. M., & Garcia-Berro, E. 

2009, Journal of Physics Conference Series, 172, 012075 
Corsico, A. H., Althaus, L. G., Miller Bertolami, M. M., 2006, A&A, 458, 259 
Corsico, A. H., Althaus, L. G., 2006, A&A, 454, 863 

Corsico, A. H., Garcfa-Berro, E., Althaus, L. G., Isern, J., 2004, A&A, 427, 923 
Corsico, A. H., Althaus, L. G., Benvenuto, O. G., Serenelli, A. M., 2002a, A&A, 
387.531 

Corsico, A. H., Benvenuto, O. G., Althaus, L. G.. Serenelli, A. M. 2002b, 

MNRAS, 332, 392 
Corsico, A. H., Benvenuto, O. G., 2002, ApSS, 279, 281 

Cox, J. P. 1980, Theory of stellar pulsations (Princeton: Princeton University 
Press) 

Driebe T., Schonbemer D., Blocker T., & Herwig R, 1998, A&A, 339, 123 

Dziembowski, W. 1971, Acta Astron., 21, 289 

Fontaine, G., & Brassard, P. 2008, PASP, 120, 1043 

Garcia-Berro, E., Torres, S., Althaus, L. G., et al. 2010, Nature, 465, 194 

Haft, M., Raffelt, G. & Weiss, A., 1994, ApJ, 425, 222 

Hansen, B. M. S., 2005, ApJ, 635, 522 

Hermes, J. J., Montgomery, M. H., Winget, D. E., et al. 2012, ApJL, 750, L28 
Iben, I. Jr.. & Tutukov, A. V., 1986, ApJ, 311, 742 
Iglesias, C. A., & Rogers, F. J., 1996, ApJ, 464, 943 

Itoh, N, Hayashi, H., Nishikawa, A., & Kohyama, Y, 1996, ApJS, 102, 41 
Kalirai, J. S., Bergeron, P., Hansen, B. M. S., et al. 2007, ApJ, 671, 748 
Kawaler, S. D., & Bradley, P. A. 1994, ApJ, 427, 415 
Kawka, A., & Vennes, S. 2009, A&A, 506, L25 



13 



A. H. Corsico, A. D. Romero, L. G. Althaus and J. J. Hermes: Pulsating low-mass white dwarfs 

Kepler, S. O., Kleinman, S. J., Nitta, A., Koester, D., Castanheira, B. G., 
Giovannini, O., Costa, A. F. M., & Althaus, L. G, 2007, MNRAS, 375, 
1315 

Kilic, M., Brown, W. R., Allende Prieto, C, et al. 2012, ApJ, 751, 141 

Kilic, M., Brown, W. R., Allende Prieto, C, et al. 201 1, ApJ, 727, 3 

Kilic, M., Brown, W. R., Allende Prieto, C, Pinsonneault, M. H., & Kenyon, S. 

J. 2007, ApJ, 664, 1088 
Kleinman, S. J., Kepler, S. O., Koester, D., et al. 2012, ApJ, submitted 
Magni, G, & Mazzitelli, I. 1979, A&A, 72, 134 
Marsh, T. R., Dhillon, V. S., & Duck, S. R. 1995, MNRAS, 275, 828 
Panei, J. A., Althaus, L. G., Chen X., & Han Z., 2007, MNRAS, 382, 779 
Romero, A. D., Corsico, A. H., Althaus, L. G, et al. 2012, MNRAS, 420, 1462 
Romero, A. D. 2012, Ph.D. Thesis, University of Buenos Aires 
Sarna M., Ergma E., & Gerskevits-Antipova J., 2000, MNRAS, 316, 84 
Siess, L. 2007, A&A, 476, 893 

Silvotti, R., 0stensen, R. H., Bloemen, S., et al. 2012, MNRAS, in press 

Steinfadt, J. D. R., Bildsten, L., Kaplan, D. L., et al. 2012, PASP, 124, 1 

Steinfadt, J. D. 2011, Ph.D. Thesis, University of California 

Steinfadt, J. D. R., Bildsten, L., & Arras, P. 2010, ApJ, 718, 441 (SEA10) 

Tassoul, M., Fontaine, G, & Winget, D. E. 1990, ApJS, 72, 335 

Tremblay, P.-E., Bergeron, P., & Gianninas, A. 2011, ApJ, 730, 128 

Unno, W., Osaki, Y, Ando, H., Saio, H., & Shibahashi, H. 1989, Nonradial 

Oscillations of Stars (University of Tokyo Press), 2nd. edition 
van Grootel, V., Dupret, M.-A., Fontaine, G, et al. 2012, A&A, 539, A87 
van Kerkwijk M. H., Bassa C. G, Jacoby B. A., Jonker P. G, 2005, Binary Radio 

Pulsars, ASP, Eds.: F. A. Rasio, I. H. Stairs, 328, 357 
Winget, D. E., & Kepler, S. O., 2008, ARA&A, 46, 157 
Wu, Y., & Godreich, P. 1999, ApJ, 519, 783 



14 



